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Chapter I 



Introduction 

Problems frequently encountered in complex multiple input-ouput systems 
stem from interaction or coupling between the controlled variables. For ex- 
ample, in the lateral control of an aircraft, a change in the angle of bank by 
aileron deflec^on will result in a change in yaw angle provided that there is 
no rudder input to offset the change. To eliminate this interaction between 
the controlled outputs, the theory of non-interacting controls was developed. 

For a control system designed in accordance with this theory, a desired 
change in one output variable will not cause a change in the other output 
variables. 

The general method for designing such non-interacting control systems 
was first developed by A. S. Boksenbom and R, Hood, Ref. 1, and applied 
by them to the control of a turbo-jet engine. This work was later repeated by 
Ho S. Tsien, in Ref, 2, In Ref, 3, H, Freeman extended the theory to obtain 
a desired degree of interaction in an m input and n output system. In Ref. 4, 
Freeman discussed stability and realizability but he did not discuss the practical 
design problems or use the more realistic open loop approach. 

This study will develop the basic criteria for non-interaction of the 
variables in an n input and n output system. In addition, the effect of dis- 
turbance inputs and errors in mathematically describing the physical system 
and synthesizing the controller will be discussed. It will also introduce the 
potential designer of a non-interacting control system to some of the practical 
problems which are often encountered in applying the basic non-interaction 

i 



criteria. As a final step, the theory will be applied to the design of an auto- 
pilot system in order to illustrate the design techniques and the difficulties 
that may be encountered in a practical application. 
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Chapter II 
General Theory 

In this chapter, the basic relations between controller and controlled 
elements of a system to effect non-interaction of variables will be developed. 

To achieve this, a specific system configuration will be assumed. However, 
the extension ®f this theory to other configurations is feasible and should be 
apparent from the subsequent development. 

Introduction 

A linear physical system with n inputs, , and n outputs, Lj-^ , can be 
simply described in vector notation. In such notation, and ij are con- 
sidered as n dimensional vectors with components and . The coefficients 
of the equations relating ^ to <^ in the physical system form an n by n 
matrix, A , which is composed of elements a.j. Note that the a^ are oper- 
ators expressable as ajj(p) where p = d/dt Thus, 

C|- = i = i toTl (2-1) 

^ J=i 

or in vector-matrix notation, 

y (2-2) 

In a block diagram, the system described above can be represented as: 










‘i 



In the case where there are m inputs and n outputs and m is greater than n, 
the m-n additional Inputs are fed directly through the system to become system 
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outputs. This implies that the A matrix is now an m x m matrix with the 
added O-n elements having a value of unity. Such a case is described at 
length in Refs. 1 and 2. 

To achieve non-interaction, it is desired that x where % is an n 
dimensional input vector with components .x.. In designing a linear system 
to satisfy the above input-output requirements, the foremost problem is to 
use the error vector, 6 = (y-y, to determine ^ for non-interaction and at 
the same time ensure that G'^0. One possible system configuration for 
achieving these ends is shown in Fig. 1(a). 

In order to keep the treatment of the subject sufficiently general, servos, 
transducers and controllers have been included in the overall system shown 
in Fig. 1(a). The transducers have been supposed in light of the fact that the 
real outputs, iji , of a physical system can not be obtained directly. Thus 
the measured outputs, related to the real outputs, , through the 

transducer matrix, T. This matrix is a diagonal matrix with elements t^^^. 

It is further supposed that the transducers are selected such that the 
are very small and hence no attempt will be made to correct to Lji , 

The physical system inputs, Si, , are generally servo driven and the are 
related to the servo signal inputs, , through the servo matrix, S. Like 
the T matrix, S is a diagonal matrix with elements s„. The controller, which 
relates the system input signal vector, S » to the error signal vector, ^ , 
is represented by an n by n matrix, C. It will be designed such that the 
overall system is non-interacting and has desired response characteristics. 

As a means of simplifying the ensuing mathematical analysis, the system 
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Fig. 1 Typical Cuu Irol £yL:tem Config'ara Lions 
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shown in Fig. 1(b) shall be used. It should be noted that this is an equivalent 
system to that of Fig. 1(a) and that many multiple input-output systems can be 
put into this form. This does not mean that this is the only form possible. 
However, conditions for non-interaction ai'e more complex and system design 
becomes exceedingly difficult, if not impossible, when configurations different 
from Figs. l(af and 1(b) are used. A simple exception to the above is shown in 
Fig. 1(c). Here the diagonal matrix R, which might for example provide deriva- 
tive feedbacks, has been added. Although another loop has been added, the 
system of Fig. 1(c) reduces to that of Fig. 1(b) if A = (T + R)AS and ij = (T-f (?)^. 
Analysis of the Control System 

The system under consideration is the simplified but equivalent system 
shown in Fig. 1(b). The vector equation for this system is: 



(j= = ACe = AC (x-^) 



(2-3) 



Solving equation (2-3) for the system output vector, ij , provides: 

1-1 



u- [l+AC]' AC X 



(2-4) 



where I is the identity matrix. 

In order that the inverse matrix [r+Ac] shall exist, it is required 
that the determinant, |H-AC| , is not identically zero over all values of the 
complex variable p. Of course, AC will be a function of the operator p and 
hence |r+Ac| may be equal to zero for a few particular values of the operator 
p considered as a complex variable. 

If the system under consideration is to be non- interacting it is required that: 
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where D is a diagonal matrix with elements d-. This also implies that 

iji = 4- Xi for all i or that a specific output is affected only by a corresponding 

input. * 

From equations (2-4) and (2-5) it follows that: 



[t-hacJ'ac = D 



(2-6) 



Certainly, one condition on the matrix AC that satisfies the above equation is 
that AC is a diagonal matrix. If AC is diagonal then [r+Ac] AC is given by; 



[i+ac]’'ac = 



(gc)„ 



O 

(ac) 






O 



l+(<ac)u 

0 

« 

0 • 



o 



where (&c)_ are the components of AC. 

Solving equation (2-6) for AC yields.- 



O 



(aC)LL 

^ I4<a£)u 



(2-7) 



AC = [i-3T'd 



(2-8) 



Again, it is obvious that AC is a diagonal matrix if equation (2-8) is to be 
satisfied. 

In summarizing, a necessary and sufficient condition for non-interaction 
is that the matrix AC be diagonal; i.e. AC - D. Furthermore, the elements 



i 
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(ac) or d must be finite and not identically minus one. This assures that; 
ii ii 

Mt=- ~~ X: = /or U I ten (2-9) 

' I'KOck ,4-da 

where d^^^ is finite. 

It is now known that AC = D where D is a diagonal matrix different from 
the diagonal nratrix D. Solving for C by matrix methods yields; 



C = A''d 



( 2 - 10 ) 



For A ' to exist, the determinant of A must not be identically equal to zero. 

This requirement is also necessary for n inputs to specify n outputs independently. 
Thus, a condition where IAN O implies a lack of independence among the 
system equations and the impossibility of completely controlling n outputs 
by n inputs Si • 

Equation (2-10) can be rewritten as; 



Cii = 



a'L d 









( 2 - 11 ) 



D is a diagonal matrix with elements d^^ and therefore the terms in the summation 
are non- zero only for k = j. Thus, 





( 2 - 12 ) 



Considering the case where j = i, equation (2-11) reduces to Cj^ =• CCiidu or; 




Cu 



a-' 



u 



(2-13) 



Note that equation (2-13) is the open loop transfer function of the overall system. 



Here is known and c_ is designed so that the closed loop transfer function, 



dii 



l+dtl 



, provides favorable response characteristics. 



After the c^^ have been found by conventional compensation methods, the 
off diagonal elements of the controller matrix, Cjj, are found by combining 
equations (2-12) and (2-13) to obtain: , 



Cii= 



. ti = — .rj. „ -for j ^ ' 

ar‘- U- 



tcp/L 



(2-14) 



If any of the OCjjEO, then the choice of djj is arbitrary. It follows that 
the off diagonal components, c.j, are then determined from equation (2-12). 
Although d^.^ is arbitrary when Cijj - O , it should be chosen such that good closed 
loop response characteristics are obtained. Thus, it should be chosen fairly 
large in order that outputs approximate inputs with a minimum steady state 
error. This may be seen from the relation, • Naturally, the d^^ 

can not be overly large otherwise the c_ must assume impractical proportions 
through equation (2-12). 

In summary, the c^^ are determined by the design of the open loop transfer 
function, d|^ = — where CC n is known. Thus the c.. are essentially compen- 

(X'a 



diL 

sabLon so chosen that is stable with favorable response characteristics. 

i+du 

The off diagonal elements of the controller matrix are then found from equation 
(2-14). Provided no difficulties are encountered in compensating the d^, the 
system performance should be as desired, namely, non-interaction of control 



variables. 



m 
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Chapter III 
Error Analysis 

There are essentially two basic types of system errors. The first type 
of error is a dynamic error resulting from random disturbances entering the 
system. The second type are those errors introduced into the system through 
the A and C mStrices due to inaccuracies in describing the true physical system 
and/or inaccuracies in the synthesis of the controller elements or transfer 
functions. Errors of the second type could obviously result in some system 
interaction and conceivably stability problems. 

Disturbance Errors 

In analyzing errors of the first type, the disturbances were considered 
as input signals to the A and C matrices as shown in Fig. 2. It is supposed that 
m disturbance signals, LLc ^ , are inti’oduced into the C matrix and r disturbance 
signals, , are introduced into the A matrix. These input signals may be 

represented as an m dimensional vector, Li , and an r dimensional vector, V . 

As an example of the source of these inputs, the Li type of input might be 

2 “ 2 - 
considered as 7T noise or drift type signals introduced at the inputs of the >T 

Cjj components. The V type of signals are introduced into the^ physical system 

i 

being controlled and hence could be torque loads, gust effects, etc. Regai’dless 
of their source, it is desirable to minimize the effect of the Li and \/ dis- 
turbances. 

It will be assumed that Li and V are related to the system in vector- 



matrix form as follows-. 






I 

I 







ii 



% 




Fig. 2 Introduction of Disturbances into the 



Control Eysteni 
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^ = C6 +- C*LL 

U- AS -h A* V (3 

In the above equations, C* is an n by m matrix and A* is an n by r matrix. 
Noting that £= Xrlj , the system error equation can be found from equations 
(3-1) to be: * 

6= + [x+ACy'(><-'^C*U.-AV) 



Specifically, the errors due to the LL and V inputs can be expressed as: 

€u= - ' AC*U 

(3-3) 

€v - ~ [x + AC] A*v 



In designing a control system, the designer will have to investigate the 
errors resulting from the U. and V disturbances. This will require an investi- 
gation of the above vector equations to determine the error associated with the 
specific input disturbance. This may be illustrated by considering the effect 
of a disturbance input to a single component of the C matrix in a two-dimensional 
system. Assuming the disturbance as an input, U, , to the c ^2 component of 
the C matrix, equation (3-3) becomes: 



e>^.=-cx+Acr 








or 



I 

I 
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As developed previously, for non-interaction, AC is a diagonal matrix which 
was denoted as D. Therefore, 



- 



r * ^ 

< +dii 






^ i+d*2 




(fl2i Ciz) 



U., 



and 






•u. 



la 



I "id I 



U., 



e,.. = - 






0.2 1 ClZ 
I -idzz 



tt, 



For illustration purposes, Gf will be the only error considered. From 



equations (2-13) and (2-14); 
dii = 



; C,2. = Cz2 = - C 



Therefore, 



n 


OrL 


a„ 


G, = 


- C„ 




O'!, 


CLzz 




CLs. Czz 







2.Z 



'tt. 






Assuming the hypothetical system represented by: 



t 



P 



A 




pCp-H) 



I 

p+l 



J_ = JAL ^ ^ _L • 

<r,', (Uz p rpF I ) p ' 



I _ p ( p i". ‘I ~ I 
pCp + l)"' ' pd-l 



then, 



14 



Now, if the control matrix diagonal elements are designed as: 



c„= 1 






then, 



C( 2 := ~>Z 



L 



p(p+')_ 






Z 



Therefore, 



6 



p L ^ p+ 



iV 



a. 



Since ^ = then the output with assumed as zero, will be: 




It is readily apparent that this transfer function presents an undesirable 

situation, since at low frequencies the magnitude of the output will be large 

because of the term .'2/p. In addition, the output diverges with time for a 

step input, U, . Thus, in this hypothetical case the designer would be forced 

! 

to modify the controller elements in order to provide good system response 
and at the same time minimize the disturbance errors, which may be difficult, 
if not impossible. 

One method of reducing the effect of the disturbances which may prove 
effective in a practical design is to introduce a high gain prior to the disturbance 
inputs. Fig. 3(a) illustrates a basic circuit for a scalar system (one of the non- 
interacting subsystems for example) where, with X~0 , the transfer function 



relating the output, U , to the disturbance input, LU , is: 



4 
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Fig. 3(c) 



Fig. 3 A Metho'i f'or Reducing the Effect of 

Disturbances Entering the Control System 
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ii = JIM- = Y(p) 

>J- i + Y,(p) 



(3-4) 



With the introduction of a high gain element yOi. as in Fig. 3(b), the circuit 
can be reduced to tliat shown in Fig. 3(c). From this latter figure, with X=0 , 



U, i+T(pj^ 



(3-5) 



By proper design, can introduce high gain over the desired frequency 
range and consequently reduce the output due to the disturbance inputs. yU. 
is, of course, designed to maintain good system characteristics for Lj with 
X as an input. 

Component Errors 

To analyze the second type of error it will be necessary to refer to 
equation (3-2) which is repeated below: 

t-l 



6= [x+Ac]' (y.-Ac*a-A*v) 



(3-2) 



The elements of the A, A*, C and C* matrices can be redefined in terms of a 
lowest common denominator as follows: 



II 


0^l\ (p) 


da (p) 


II 


dc) (p) 


^ (p) 


C(j = 


'a. 

‘M) 


dc (p) 


It 


^ Jk 






(3-6) 



I 

■ 



I 



) 
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where , Ccj , (Li.j , C-ij are polynomials in the operator p and dj^(p) and 
dg(p) are the lowest common denominators of the a^j and c^j components re- 
spectively. 

From the above definitions, it is apparent that the al have the same poles 

as the a., by definition. The same is true for c* and c... This generally is the 
ij ij 1] 

case and if it is*not, then the nonacceptable poles may be removed and used to 

! 

define a new vector, (X or V • It is now possible to define new matrices 
whose elements are polynomials in the operator p. Thus, 



X-(a.j) ; A‘=(a-) 

(3-7) 

C.(q) ; C- = (CJj) 

Note that /\ - A* = ~ f\* ^ C = ^ Hence, 

OcL ^ da dc dt 

equation (3-2) may be written as: 



£ = l^do-cjc I -f-AcJ X ^ ~ A C U,-dtA v) ^3_3J 

In the above equation, the components of A , A , C , C , and dc 

are polynomials in the operator p. Hence any poles of equation (3-8) result 

I I -T— ^ 1 I i ^ I 

0(j.ac X A C . If determinant Oadc X 4- A C j 

has no zeroes in the right half p plane, or on the jw axis (imaginary axis), then 

the system is stable for disturbance inputs as well as normal inputs, Xi • 

i \ I 

ao.<)cX-l"AC and d d have zeroes which are 
' a c 

common, and if these zeroes are in the right half p plane then the system will 
be stable for X inputs but could be unstable for (X and V inputs. If d^d^ 
has no zeroes in the right half p plane and the system is designed so that 
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diz ^ 

, .1 and 777- are stable, then the error vector, cz , must be stable for all 
' '*'^11 l + Oii 

inputs, X , U. and V . 

When there are small changes in the coefficients of A and C due to inaccur- 
acies in the synthesis of C or inexact knowledge of the coefficients of the equations 
of the physical system, the zeroes of the determinant of -t-AC will 

change because t)f the coefficient change. Hence, the poles of equation ( 3 - 8 ) 
will be altered slightly for small changes in the coefficients. Thus, if the 
system is designed stable and if d„d has no zeroes in the right half complex 
plane we are assured that the system remains stable. 

The small changes in the coefficients of A and C will result in interaction. 
However, if the deviations are truly small, then the resulting interaction will 
be negligible. This is due to the fact that the d^j or off diagonal elements will 
be small and hence will have little effect on the overall system. Thus, it may 
be seen that interaction due to these slight deviations can be minimized only 
through careful attention to the analytical description of the physical system 
and accurate design and mechamzation of the controller. 



! 



> 
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Chapter IV 

Practical System Problems 

In the previous chapter, the effect of small deviations, ^ , in a.j 
and c^j were discussed. In this chapter, the problems or consequences 
attendant with gross deviations in c^j because of design advantages will be 
investigated. Ig addition, problems arising from impossible or impractical 
compensation requirements, as may be dictated by non-interaction theory, will 
be discussed and means of alleviating the problem pointed out. 

Problems arising from gross deviations in c -j. 

Suppose that in synthesizing complicated controller transfer functions, 
c^j, simplifications are made by neglecting terms in the numerator and 
denominator of c^j that have "break" frequencies well above the natural 
frequency of the system. Further suppose that all these discarded terms have 
negative roots. The obvious consequence of such gross deviations from the 
desired c- is system interaction at the higher frequencies. In addition, one 
must consider the effect of these departures on the stability of the overall 
system. 

To investigate these ideas more fully, assume a two-dimensional control 

system in which the exact controller transfer functions, Cj^£ ^Jid C21 are 

altered by a factor q to produce new relations c’ and c' . Analytically, 

12 21 



c 

c 



/ 

/z 



= f' 



C 



/Z 



(4-1) 






(4-2) 



I 
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For the theoretical non-interacting system, the following equations have been 
shown to exist; 

- ^11 ^'2 Czi = X| = (4-3) 

^'2 ~f~ 0^22. (^22 ^ (4-4) 

Cz -h diz Ciz = 0 (4-5) 

^2/ = ^ 2 i Cl, -i- dzz Czi =: 0 (4-6) 

Substituting equations (4-1) and (4-2) into equations (4-3) through (4-6) where 
applicable and manipulating and rearranging where desired, the results are: 

dll = dll Cl -i- Ci,Z (^2. Czi = 

d^ - duCzz 4 - (Iz! Ciz Cz, Ciz^JCzz (4-8) 

d,z ~ G-n Ciz -hdix Czz ^ j^d/z 0 ( 4 - 9 ) 

dzi — 0.ZI Cn -h (^2 Czi dzz ~ j^Clii Cj ^ 0 (4-10) 

where the d' are elements of a new and non-diagonal matrix, D'. 

ij 

The interaction is evident from equations (4-9) and (4-10). The changes 
in the transfer functions d^ due to the factors q^ may be determined from 
equations (4-7) and (4-8). 

In investigating the effects of gross deviations on system stability, it is 
best to look at the system equation m vector-matrix notation. It is; 



i 






t 
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y = [l -hD'l d' X 



where D' is a non-diagonal matrix. 

The stability of the closed loop system is determined from the determinant 
of jXi-D I which the characteristic equation of the overall system. It can be 
shown that: 



ll+D'l = (n-d'„)0+d^)-d,id':i 



(4-11) 



It follows that if this resulting polynomial has zeroes in the right half complex 
plane the overall system is closed loop unstable. 

One means of determining the stabiliiy of the closed loop system without 
solving for the roots of the polynomial is as follows: 

( / + dll ) ( I fdzz ) - d,2 dzi = 0 



or when rewritten: 



(l+d;,)(i+d22) 



r 



dU, 



(4-12) 



stability' is assured if none of the above factors have zeroes in the right 
half p plane. 

The open loop transfer functions, d^^ and d^^, which are given in 
equations (4-7) and (4-8), should be plotted in Nyquist plots and Nyquist stability 
criteria, N = P - Z applied to determine if there are zeroes of the 1 "h d|j In 



the right half p plane. 



I 
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Assuming that there are no zeroes of the 1 + d! . in the right half p 



plane and that d '^2 d^j^ have no poles in the right half p plane, then for 
system stability, the absolute value of , — -rrrr -r — r must be less than 



(l+duX • +diz ) 

unity for all frequencies, tO . However, reference to the Nyquist plot of 
a typical d|^ shown in Fig, 4 provides the following important relationship; 






Ll 



a 



where 



a- 



Mi-hl 



(4-13) 



It follows from equation (4-13) that: 

1/ i' 

( ^ 



ctiz dll 



dl'zdi, 



(i+M,)()+f^0 (l ^ciii )(l+dzX) 



Thus, if (1 + M^) (1 + Mg) d^g d^^^ < 1 for all frequencies, LO , and (1 +- d^^^) 
as well as (1 +• dgg) have no zeroes in the right half p plane, then the overall 
control system is closed loop stable. 

Compensation Problems 

It is conceivable that the co-factors of the a^ of the physical system 

may have zeroes in the right half p plane so that the ~~~r will have poles in 

CLii 

the right half p plane. From the basic open loop non-interaction relation, 

, it follows that it may be impossible or impractical to compensate 

is stable. In this event, it would be impossible 



I 

OTi'l 



d“ 

with c.. so that — !i- 
IFdu 



to obtain a practical non- interacting control system using the design techniques 
discussed in the theory section. As an example, consider the following simple 
system: 



A = 



1 


. 2 , 


P+1 


p+l 




-.Olpt 1 




p^l 



: ^-(V 
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Fig. 4 Nyquist Plot of a Typical Open Loop Transfer 
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(~.0I p-hl)(p-hl) 
01 p 9C> 



- -• 01 P -f-' 

i-.oip ^/)(pi-l) 

The locus of — ^ for p z jw is: 

^-1 




But for stability, N = P - Z and in tliis case P = 1. Therefore, N must be 
if Z is to be zero. However, from the general non-interaction theory, d,,= 

It is quite obvious that the compensation of such that N = -I- 1 and Z = 0 is 

a difficult problem. One means of compensation is to use the prediction oper- 

p 

ator, € . This will provide the necessary encirclements but is entirely 

impractical as a means of compensation in a control problem. Other than this, 
no compensation could be found. One might consider mechanizing the compensa- 
tion so that it cancelled the objectionable poles and zeroes. However, this wiU 
not work because the poles and zeroes of the physical system are not known 
exactly and unless we have exact cancellation the poles and zeroes are not 
really cancelled. Thus, it appears that the straight forward design techniques 
develo|)ed in the theory can not be used when confronted with such a compen- 



f 1 

a-.! ” 



Then; ^ ^ 

(A, II 



and; 



a- 



- 1 
II 



sation problem. 
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In some cases the foregoing design problem can be avoided by merely 
interchanging the rows of the S vector. This implies an interchange of 
controls. If the design problem is removed by such a change, one may continue 
to design the control system in accordance with the general non-interaction 
theory. 

As an exaAple of the foregoing, consider the same system of the previous 
example but with the ^ vector rows interchanged to form a vector S and a 
corresponding change of the columns of A to form a new matrix A'. Thus, 




-•Olp-H 



_ p-H p-H 

The new d- are now: 




C.KI = 



- Cm 

P4-I 



^22 = Yo. = ^^1^1 ~ CZZ. CZ) (- Olp +-,9(p) 

^ / I 



There are no poles in the right half p plane of the above open loop transfer 
functions. Hence, P = 0, In addition, the design of c^^j^and C 22 so as to have 
N = Z = 0 in the closed loop system is not anticipated to present any real 
compensation difficulties. Thus, the impossible or Impractical compensation 
problem has been avoided. It should be emphasized that interchanging the rows 
of the ^ vector will not alleviate this particular compensation problem in all 



i 








I 

I 
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instances. This is obvious from the foregoing example if either Tty or 

0.|2 

I 

Tjrr have one or more poles in the right half p plane. Naturally, the presence 

of such poles does not imply that the can not be practically compensated 

o.~q 

so that Z = 0 but rather that difficulties in compensation are more probable 



than the instances when P - 0. 



In a more Complicated system where there are n control inputs and out- 
puts and n is greater than two, the poles of — i — result from the co-factors of 

^ ll 

the A matrix. These complicated co-factors tend to increase the probability 
I 



that wiU have one or more poles in the right half p plane. Thus, one 

^ d 

might expect some difficult and impossible compensation problems when applying 
the general non-interaction theory to complicated physical systems. If an inter- 
change of the controls, (interchange ) does not avoid these problems then a 
different design approach must be taken. 

The theory behind one approach and the effect on system performance is 
most easily described for a simple two-dimensional system. For instance, 
suppose that d^^j^ can not be properly compensated by The result is that 

the system can not be designed entirely within the non-interaction theory. How- 
ever, a practical system can be found if the designer will accept some inter- 
action. Furthermore, some of the non-interaction techniques can be used once 
the difficulty has been removed. 

In finding a practical system it is probably best to start with the basic 
expression for This is: 



du — dt\ C(| •1' 0-izCzi 



(4-14) 
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In as much as it is impossible or impractical to compensate with c^j^, 

then Cj^j^ may in general be chosen arbitrarily. However, in some instances Cg^ 

may be chosen completely arbitrary and hence the arbitrariness of Cj^j^ will be 

restrained to a certain type of transfer function which results in 1 + dlj^^ having 

no zeroes in the right half p plane. In addition, certain simplifying assumptions 
I 

in may re^lt in restrictions on In any event, c^^j^ assumes a new 

description as Since it may be necessary to change C 2 ^ from the transfer 

function dictated by non-interaction theory, it now becomes c^^- Therefore, 

I 4-d|( = I ^11 Cn 4- 0.,2 Cu (4-15) 



The procedure is to make 1 f so that there are no zeroes in the right 
half p plane. Note that 1 -f- dj^^^ will have zeroes in the right half p plane because 
dfi = can not be properly compensated so that Z = 0. 

In as much as the a^ are fixed by the physical system it is much easier 
to design c' and/or c' so that there are no right half p plane zeroes in the 

JL JL ^ X. 

resulting 1+ d.' . As a first approximation, write c' in the form dictated by 

-i- <l» ^ X 

non-interaction theory, which is: 



Cai — 



<Xzi 



CLzz 



(4-16) 



In as much as c^^^ is arbitrary in some sense, make it as simple as desired. 

Then substitute into equation (4-15) and apply Routh's criteria to the resulting 
polynomial to determine whether or not these are zeroes in the right half p plane. 
If this test fails, alter c^^^ by a factor, a , where a is a polynomial in the operator 
p or a ratio of polynomials. Again, test the resulting polynomial for positive 
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roots. Eventually, a q will be chosen such that 1 + has no zeroes in the 



right half p plane. In addition, if q is chosen so that the "break" frequencies are 
above the system natural frequencies then the interaction will be restricted to 
these higher frequencies. ! 

Having removed the difficulty in d , attention can be focused on d„„. An 
attempt should b^made to design it in accordance with the general theory of non- 
interacting systems. If no compensation difficulties are experienced then c^g 
and Cj, 2 will be dictated by the general theory. In addition, 1 d 22 = 1 +• <^2 



and there will be no right half p plane zeroes. If compensation difficulties are 
encountered then follow the procedure outlined for d!|^j^. 

Naturally, the system has some interaction which appears through the 
d !|^2 *^21 functions. In fact, for the case presented here 



Due to interaction, the overall system closed loop stability is dictated 
by the zeroes of: 



Determining the stability through this equation has already been discussed at 
length in the previous section. In the event that there were zeroes of equation 
(4-19) in the right half p plane they would stem from the interaction terms, 
d'j ^2 3Jid d^j, . Thus, it would be necessary to choose another q factor so that 




(4-17) 



dj-j = ^2(C/ 4- (X. 12 . ^2.t ^ ^ 



(4-18) 




(4-19) 
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'I A 

both the 1 -f- d!. and contributed no right half p plane zeroes to 

equation (4-19). 

In the case considered here diz— O . Thus, sj'Stem stability is assured 
if ( I 4- d(( ) ( I +■ dzj. ) have no zeroes in the right half p plane. Actually, 

dll = 0 due to the slight deviations in a^^ and of equation (4-17). How- 
ever, it is not expected that such slight interaction through this term would cause 
instability when coupled with d^^. 

In summary, when compensation difficulties arise and interchanging the 
rows of the S vector does not alleviate the problem, then it is necessary to 
deviate from the non-interaction theory to the extent that there are no zeroes 
in the right half p plane of the resulting characteristic equation of the system, 
namely: 



+d|i )0'^dz2.) — d/zdz* —o 



Naturally, interaction will be introduced through the d| . transfer functions. In 
the case considered here, the amount of interaction as a function of frequency 
can be found from the following transfer functions: 

_ d;. 









‘/.A = ik 

/X, p,| 

[JDI = (i + dii) — di, 



(4-20) 

(4-21) 



where 



I 
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Chapter V 

The Synthesis of an Aircraft Control System 
In present day aircraft, the auto pilot performs many functions in the 
integrated control system. Although the system inputs are generally restricted 
to combinations of elevator, rudder, aileron and throttle, the system controlled 
outputs are many? For instance, an auto pilot system could control the aircraft 
velocity vector through the control of the pitch, yaw and roll rates. Other 
systems could control either altitude, angle of attack, airspeed, glide path, 
Mach number, lateral acceleration, turning rate, etc. , or combinations of the 
foregoing. It is the intent here only to point out that the auto pilot can be used 
to control many variables and each subsystem of the auto pilot can control as 
many outputs as it has inputs. 

In order to point up the theoretical analysis and the practical problems 
associated with a design based on the general theory of non-interacting controls, 
an auto pilot, subsystem shall be synthesized in this section. Naturally, the 
goal of such a study is to arrive at a controller design that permits good system 
performance as well as the desired non-interaction. However, as was pointed 
out in the preceding sections, a perfect non-interaction over all frequencies is 
not obtainable nor is it necessarily required. In addition, complex physical 
systems such as an aircraft could easily require impractical or impossible 
controller compensation when applying the general non-interaction theory. 

Thus, the designer may have to accept some interaction in order to preserve 
system stability. Furthermore, he may be forced, at least in part, to deviate 
from the general theory in order to design a practically realizable system with 
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desirable characteristics. 

The subsystem to be synthesized is one dealing with longitudinal control 
in which the inputs are elevator and throttle deflections. The controlled out- 



emphasized that this study was chosen only as a means of lending practical 
significance to the theory and is not considered as a demonstration of a specific 
auto pilot design. 

The assumed aircraft and flight condition is an F-86 "Saberjef making 
an approach down a 3° glide slope. Small perturbations from symmetric 
equilibrium flight will be assumed thus permitting the use of the linearized 
longitudinal equations in the ensuing analysis. Upon determination of the 
controller design for the best system performance, the synthesized system 
will be studied on the differential analyzer. Such a study should clearly demon- 
strates the applicability of the theory, error analysis and the practical system 
problems. 

The large motion longitudinal equations of motion for an aircraft using 
a combined body- wind axis system are shown in Ref. 5 to be; 



puts are glide path angle, , and indicated airspeed, Afp . It should be 




(5-1) 



-h XI5 — yy)^sin&^ 

^ Rwa — 0 ~ ~)ToV’p(^ I'Rvi' ~ — Fx -k COSoi 



( 5 - 2 ) 




=0 Q +[r»-ixjRP 

+ -t- rx^ (p^-QR) -i-Ej? ('r-pq') 



(5-3) 
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The linearization of these equations is carried out in Appendix I. In addition, 
a fourth linearized equation is obtained so that there are four simultaneous 
equations involving the four output variables, AJp , the inputs, 

Sj- and Sq . 

It is desired to determine the transfer functions that relate Alp and 
Sj- and , It stiould be clear from the theory that these transfer functions 
are the a^ elements of the K (aircraft) matrix. The mathematical steps taken 
to determine these elements from the four simultaneous equations are completely 
outlined in Appendix I. 

The elements, a^j, are merely a ratio of polynomials in the operator p 
or D. However, they derive their existence from the coefficients of the line- 
arized longitudinal equations. These coefficients are composed of thrust terms, 
gravity terms and aerodynamic force and moment terms expressed through the 
aircraft stability derivatives and trim lift and drag coefficients. 

The stability derivatives for the F-86 were obtained from Ref. 6 for 

flight at sea level at 163 knots. They are listed in Appendix II along with 

other physical data pert aining to the F-86. Calculation of the trim lift and 

drag coefficients as well as the b.. and f coefficients of equations (1-9) 

ij ij 

through (1-12) are also found in Appendix II. Knowing the b- and f^j co- 
efficients permitted the calculation of the a- elements from equations (1-18) 
through (1-21) of Appendix I. The results are: 

-h3A77l) 

(B^F|.4mD l-3.45'34)(DV.0i5<bD f. 024 S) 



a, I = lO.ZZ 
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. 885“ 



(I>-iSo)C3>-H.|) 



CL2,\ = . 00(o3.5 



CioJL = .07 



^J^\A mj> 3.4534 ) ( dV. oisfc D +■• 0^4^ 

fP-l- l.40q) (D'"-K70ID -1-3 .087) 

DV(.4I44P 4- 3 . 4534 ) (dV.O/58D +.0245) 

(J)-. 007J7)(l^-f'7.i>5'){l>-8.S) 



(' D^4- 1 .4 144 P -I- 5 -4534) ( dV. 01 56X) f . 0245) 



The two quadratic factors in the denominator of the “5^ terms provide 
information as to the period and damping of the short period and phugoid motions. 
The period and damping factors obtained are: 

Short period motion Phugoid motion 



X = S2. = 



COric 



40 sec. 



'1^ = = 3. 38 sec. 

UJyis 

IS .382 Jp = .0842 

It may be seen from Figs. 1 and 5 that in vector-matrix notations: 
=. TAS S where 1” ~ ("tu) 

S =• (Sii,) 

A - (Zij) 






For the specific problem under study: 






Lj^=^ Afp ) ^ ~ J 

Cjz - 0VW ^ > ^2 ~ =■ Si 



T 

& 



Therefore: 



= (-^(1 dll 3,1 ) St diz S 22 ) Se — d,i Sy f (X,z Se 
©W — ( ^ 2 Z <^21 5 ||) Sy '^(Szz dzz 62 )) Se == ^21 Sy -h dzz.Se 



(5-4) 

(5-5) 



To determine the a^ terms of the above equations, it is necessai*y to specify 
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the elements of the diagonal servo and transducer matrices, S and T respectively. 
This has been done in Appendix III. In addition, reasons governing the specifi- 
cation of these elements or transfer functions are included in the appendix. With 
this knowledge, the a^^ terms of equations (5-4) and (5-5) are found to be; 



= 10.22 



^'2 dd>5 



(D-.0l2)(p^-H.42aD -1-5.4771) 



*— • . 



.07 



4D*-H|.4I44D -f-3.4534)CP^KC>l5fcD f.OZ4£r) ( 2 D+I) (. ZD-hl) 

(l>-)Sb) (D-H.l) 3 

^"^+1.41443) +- 3 . 4554 ) (oVoiSbD -h.O245)('.04D+-l)(-2-I>+l)^ 

(Df l.4(y))CD^-k>~70IJ> + 3.oe>7) 



_(P^-H.4I44D+3.4S34)( pV-OI 5bD f . 0 24^) f 2; D-l-l ) 

(D~.007I7)(P+-7.^5) (^D-e>.5) 

1.4144 X) f 3.4534)(DV.OI5feD4-.024S) ( • 04 D+l) 



It has been pointed out that d. . = c /a7.^ is an open loop ti'ansfer function 

11 ii 11 

for one "channel" of the auto pilot system. From Fig. 5, it may be seen that 
the open loop vector-matrix equation is Lj - AC G . It follows that; 

d,. = \ = y'/e, = ‘-"/a',', = ^ (5-6) 

& 

(5-7) 



€Ajfi 



It is evident that equations (5-6) and (5-7) are the open loop transfer functions 
for the velocity and pitch angle "channels". It is desired that these channels 
be closed loop stable as well as having desirable response characteristics. An 
attempt is made to design the controller transfer functions c and c so that 
these requirements are met. If they are satisfied, then the other controller 
transfer functions, Cj^g ^21 determined from equation (2-14). 
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The Open Loop Control System Configuration 
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The determinant of the A matrix was found to be; 



|Al 

|Al = .71 



= 1 ^ 11^22 ^ 12^21 



(J>-t-7.74KD-8.S87yD^-H.4lSDf3.4Si;>)(D\-Qia7J)-KQai^) 
(2Df I }(. ZPf l)(.04J>f I) (dV|. 4I44D 4i.4S34-f(D%.0l5fcI> 4.02.45^ 



Substituting into equations (5-6) and (5-7) and cancelling equal or 
essentially equal terms in numerator and denominator (excepting terms in 
right half plane ) provided: 



X , = Cl Oo . z ) 



( p-&.se>7 ) 



(p-.007<7)(p-8.S)(zpf|)C2p4l) 



'Joi = C2^ (1333) 



(p-.OI2)(. 298 p*’4-. 407p -HI )(.04pf I ) 

dn 



(5-8) 



(5-9) 



In deciding what c, , should be in order that — be closed loop stable 

with good response characteristics, Nyquist plots had to be used. A sketch 

Y 

of the Nyquist plot of °i is shown in Fig« 6. There are two poles in the right 

C|| K 

V 

half p plane of _12 l . Nyquist stability criteria indicates that N z +• 2 for 
QiK 

Z = 0. This implies that Cj^^^ must be chosen such that the Nyquist plot of the 
resulting or dj^^ i^ust encircle the critical point, -1, 0, twice in a counter 
clockwise direction. Without doubt, the required compensation presents a 
rather formidable problem. In fact no practical compensation could be found. 
Thus, it is clear that the system will be closed loop unstable in that it can not 
be properly compensated by Cj^^" 

In an attempt to alleviate the above design problem, the rows of the 
vector were interchanged. As a result, it was possible to compensate d^^^ in 
accordance with the general non-interaction theory. However, compensation 
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Fig. t) A Sketch of the Nyquist Plot for the Transfer 



Function, 



«./C„ K 
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found. Thus, this method failed. 

The next approach was to alter d^^^^ with a factor q such that there were no 
zeroes of 1+ d^^^ in the right half p plane. This approach is discussed at length 
in the section on system problems. From that section, it is known that a more 



It was pointed out in Chapter IV that c^j^ or ^21 could be somewhat 
arbitrary when could not be properly compensated by the non-interaction 
design methods. In the case at hand, instead of making completely arbi- 
trary, let us assume that the p-8. 5 pole is cancelled by the p-8, 587 zero in 
equation (5-8). Of course, it is understood that this is incorrect in that one 
of the right half p plane poles are being neglected. Nevertheless, compensating 
the remaining Y. with P * 1 results in c* being: 



Naturally, this is not the necessary compensation for closed loop stability of 
— as dictated by non-interaction theory. 

As a first approximation to let it be dictated by the equation: 




here: 



=■ Cl +■ UizCz, -hi 



(4-15) 




(5-10) 




Substituting for & 2 i, and c^^^ results in; 
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—.0 0^9^(^7np^-l)(.04p^l)f.3^4f>^ -h.izv^-h l) 

0 ip-f-0 ( p-'Oojn) 0/309 pi-OCni&p- 1) 



(5-11) 



Substituting equations (5-10) and (5-11) as well as aj^^ ^2 equation 

(4-15) and applying Routh's stability criteria to the resulting ninth order poly- 
ps 

nomial indicates that there are zeroes in the right half p plane. However, by 
altering c^^^ above by a factor q 2 » where q 2 = -(. 1178p - l)j{. 04p-l- 1), the 
new C 2 j^ becomes: 



(.lllptl)Q324p^+.ZZlpi-l) 

~ (p-.oo^n)O ! p-/-i) (. 130 ^ p-f-l) 



(5-12) 



Going through the same substitution procedure and applying Roth's criteria 
indicates that there are no zeroes of 1 -f-d^j^ in the right half p plane. It follows 
that the velocity "channel,” is now closed loop stable. In addition, inter- 
action is present in the system because the c'^^^ and of equations (5-10) and 
(5-12) are different from the c^j^ and C 2 j^ that would have been dictated by 
following the non-interaction theory. 

The d 22 transfer function was determined from the non-interaction theory. 
It is shown in equation (5-9), A Nyquist plot of -Y^ shown in Fig. 7. 



I.2P-H 



d 22 (or Yq ) was properly compensated so that 



By making C -t5 
22 

there were no zeroes in the right half p plane. Furthermore, the compensation 
was chosen so as to give good transient response. 

Having determined c in accordance with the non-interaction theory, then 



it was permissible to find c ^^2 from equation (2-14). Therefore: 




3 10 = 
50 ' 



300 '' 

60 



230 ^ 

^30o 



220 ® 
\ 40 '' 



1 TO'* 
1 90 ’ 



1 80 ® 

1 80 ® 



2 10 ® 
1 50 ® 



1 90 ® 

! 70 ® 



200 ® 
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C,Z = = -Z.O‘tS(lZp1-l)(-<ilpi-l)(zpM}(-OOt7i:,-l) 

" (p-.OI2)(.o4pi-l)(.Oip-hi)(.2i^p'^-i-.‘^<ip + i) 

In as much as terms such as (. 00667p - 1), (. 02p-»- 1) and (. 04p + tl) are 
only effective at the higher frequencies, c can be simplified by neglecting 
them. Naturally this will give rise to interaction at higher frequencies. But 
this is of no real consequence provided that these gross deviations do not im- 
pair the system stability through the d' d’ term of the system characteristic 

12 21 

equation. In the problem considered here the stability is not impaired by the 
interaction terms. 

In light of the foregoing: 





c^,04S 

(p-.O/Z) C 29?p'^-h. 409 p 4-1 ) 



where 




(. 04p-hl) (oZp-hi) 
- Cooh7p- /) 



(5-14) 



Attenuation plots of c’^^, c^^ and the c.. of equations (5-10) through (5-14) 
are shown in Fig. 8. One may readily see that the effect of neglecting the high 
frequency terms is to level off the gain at high frequencies, such frequencies 
being above the natural frequency of the overall system. 

The amount of interaction contributed through dl^^^ ^d d^^^ can be found 
from equations (4-20) and (4-21), In the case considered here, dl^^^ ^d d^^ 
are as follows; 




0,12. Czz. 



(5-15) 



42 







qp ^ JO/? J /A/5 V/V 



I 




43 




(5-16) 



where c' is given by equation (5-12) and c q is given by equation (5-14). 

^ -i- X ^ X 

Computer Simulation 

With the c.. designed, the auto pilot system was mechanized for study on 
the Reeves Electronic Analog Computer (REAC). The mechanization used is 
shown in Fig. 9. 



System response to step inputs of command velocity, , and command 

\J 

glide path angle, , are shown in Figs. 10 and 11. Note that the interaction 
is essentially zero except at high frequency. However, this was expected in 



that high frequency terms in c'^^ and c^^ were neglected. There Is, although 
not evident in Figs. 10 and 11, a low amplitude damped oscillation that affects 
the steady state chai'acter of and 0^ . These oscillations are the direct 
result of disturbance inputs, , that enter the system through the c„. In 
fact. Figs. 12(a) and 12(b) show these oscillations in il^and 0^ when there 

3 



are no inputs other than disturbances in tlie system. 

As a means of verifying the magnitude, frequency and damping of the low 

amplitude oscillation, a step input, X 2 , was introduced Into C 22 through the in- 

put to amplifier R-8 of Fig. 9. The effect of the disturbance input on — ^ — and 

J 

Svv is shown in Figs. 13(a) and 13(b), 

An analytical study of the system was made in order to determine the 
relation between and x and also to verify the traces seen in Fig. 13 as 
those that stem from disturbance inputs. In the case where there are no desired 
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Fig. 10\a) Aircraft Velocity 
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Fig, ll(a ) Glide path angle,© 
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Fig. 11 Syetem Response to a 10 Volt Step Input 
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Fig. 12 System Respoo&e Due to Actual Disturbance Inputs in the 



Fig, 13(a) Glide rath Angle 
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Fig. 13 System Response to a 5 Volt Step Disturbance Input to C21 

(A)p^ = =0) 
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perturbations from steady state condition, then Therefore, it may 

be shown that; 



Ajp 



(Ziz. Czi 

I -f- diiCii 'hdizCzi 






(5-17) 



Substituting into equation (5-17) and simplifying results in; 







(pi-l.l) rp-/go)(. 3^4p^-/-. 227p-hl) 

+- 44.1 p® f’^OOp'' + 5^70 p** f 2T0SDp^-»-40,4OOp*tl3lfe0Op^ + iVIfOaOp"^ 

■hSiSOp 



(5-18) 



where the static gain is unity. 

Application of Routh's criteria to the ninth order polynomial in the denomi- 
nator of equation (5-18) indicates that there are no positive roots. Therefore, 

the system is stable for the disturbance input, x . In addition, it may be shown 
2 

that (p . 0194p +• . 0210) is a factor of the ninth order polynomial. Further- 
more, it is the only factor which has a low characteristic frequency. In fact 
this frequency is . 0231 cycles per second and the damping factor is ^ = .067 . 

It is quite apparent that this quadratic factor alone accounts for the lightly 
damped high period oscillation seen in Fig. 13. Naturally, similar analysis 
can be made for other disturbance inputs into the system. 

The low amplitude oscillation discussed above as well as other responses 
to disturbance inputs can have a detrimental effect on the system due primarily 
to their magnitude and damping. The addition of a high gain loop would introduce 
gain prior to the introduction of noise and hence make the system less sensitive 



to the disturbance inputs. 
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Chapter VI 
Conclusions 

1. In theory, control systems based on non-interaction of control variables 
are feasible but somewhat complex. For complete non-interaction, the theory 
requires precise mathematical knowledge of the behavior of the system to be 
controlled aiwl exact mechanization of the controller transfer functions. In 
reality, these requirements generally can not be met but can be approached 

to the point that the resulting interaction is negligible for all frequencies. 

2. Small deviations, , hi a- due to insufficient knowledge of the physical 
system and/or in c^^ due to the inability to mechanize transfer functions exactly 
should not impair the stability of the system if the d_, or open loop transfer 
functions, are designed with a reasonable amount of stability. 

3. Gross deviations in c , so as to obtain a simpler mechanization can 

ij 

result in considerable interaction and possible instability in the closed loop 
system. If deviations are restricted to terms effective only at frequencies 
well above the natural frequency of the system, then the interaction wiU be 
restricted to these frequencies. However, system stability can be impaired 
through the interaction terms now present in the characteristic equation of the 
overall system. This equation should be investigated for zeroes in the right 
half complex plane. 

4. System errors due to disturbances entering the system can be large. 
However, these large errors can be reduced considerably by addition of a 
suitable high gain loop. 

5. The system error equation indicates that the overall system can be 
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stable for command inputs, x., and yet be unstable for the disturbance inputs, 
Uo^and . In any event, a complete error analysis for disturbance inputs 
should be made. 

6. Zeroes in the right half p plane of the co-factors of can create compen- 
sation problems through the open loop relation, = ^ , Naturally, c.. must be 

0^ U 

chosen so that there are no zeroes in the right half p plane of 1 + d... Making 

11 

the proper choice of c.. can be a formidable problem in some instances with 

11 

the result that no practical compensation can be found. 

7. Interchanging tlie rows of the control vector, S , will alleviate the 

impractical compensation problems in some cases. ‘ 

i 

8. When compensation problems are encountered and changing the rows of 
the $ vector does not alleviate the problem, the non-interaction theory can 
not be used for the design of that particular "channel" or open loop part of the 
system which is causing the problem. Some other approach must be used to 
make the "channel" closed loop stable. Normally, the end result of deviation 
from the non-interaction theoi’y is system interaction, which, it is hoped, can 



be made small. 
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APPENDIX I 



^OF raE°SajENTS 

A. Llnsarlzsd Longitudinal Equations 

Th, aquations for longitudinal notion in th, ganaral 
or large nation case in a conbinad body-^ind axis systan 
nave been shown in Rof. 5" to bo: 



(z-z) 



+ - mcvjsio^^ 

sJ 



R>/r-0-4.mVp(5^+F^g ^ 

-i- Cas&^ ■ 

=0 - -lyy Q + -Txx] R.P 4- r?x 0^ ~ p 3 4- /M 4-r 

+ rx.v Cp-hqr] ^ 

considering snail perturbations Frol a steady state 

ongitudinal Flight condition, the above equations can be 
linearized as follows: 

1. Assume the Folloving conditions: 

P=^-/3-<^=^0 * J~ -r. 

- Qw, - o(^ =z o • 7^-0 

2. Assume the Following steady Lt. plus perturbation 
relations For the system variables: 

vp=vp, 

; 4 = 4, ^ 4 

©« = av, i- ©., ; 4 = q. r 

Q - Q< 
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where Afp.P ^^A} o Q.^ Sr the perturbation 

variables.’ ' 

3. The dimensionless coefficients, Xg,Zg,Yg,Ng, Lg , and 

Mg are functions of ^ fk ^ Qs Rg« By using 

non-dimensional angular rates, Q and R , 

then Xg , Zg and Mg become functions of c><^ (X ^ 

and control inputs only. Thus, neglecting the lesser 
Important derivitives, the dimensionless coefficients 
become : 

(1-4) IXI5 = -^Cy.^ 

(X-5) 

Ms =~p^p'^ ^ +" -P Ctv)T Q 

J 

4. In the wind axis system, it is true that: 

C^-7) — Qiv 

(X-8) 0,-Qw 

5. Linearizing equations (X-l) through (x-3) about 
the steady state condition (subscript one) and 
rearranging in operator form with D * d/dt provides: 

(X-9) b), 4- b,2^ 4-b,3, 0xAy ~pU. » = X), 

vjObeni : 




= o 



IH* 



I 

I 
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Sii — 0 

J~,z — 



— CoScv'i 5. mo/, 

mjL^SrJ, ' \^i, 



(in I o) b>2_i 4Tp + ^2.2 ^'h^4 Q - 



. ‘ a. + '^^ 1 . 

t« - li- == D - p-^— Ca.v - i? COSX,-t;Sllrw’^ 



t>z 3 = 
b ^4 =i 



'/p, 

-[| 4 -K.j] :--- [M-f|£a^_ 



< 



hi ^ 

^2Z ^ 





Si n<y, 




eiii c 

_ 




) t>3i /Op -f" t>32 4" t>33 Gy 




The fourth equation is found from equation (i--7) 
and (X-S) as follows: 



— Q ~ cc which when linearized 

m 

with — ^1 — O provides: 

^ - - 

Ow ^ = o 

Rearranging and using operational form results in: 

w* - - 

(X- 12.) A^p -h b4 3 -h l>M- ^ ^42 

b^i — o 
b^-z = I) H 

^ D a d/c^t 
fc44 = —I 
- O 

~^A-Z — O 

B, Determination of the Elements of the Aircraft Matrix 
The vector-matrix equation for controlling airspeed, 
aJJ) , and the glide path angle, 0^, using power and elevator 
inputs, £j and is: 




V/hen the relation above is expanded, the result is: 

14 -) (dJp — d-n Sx -h CLii, S(>. 

(n-is) ©yj =■ Q-zx Q-ii. <Se 

The elements, CLij , in equations (X-H) and (x-f5) are 
found by solving equations (X-S) through (x-iz) simult- 



aneously for 4^ and • Using determinant methods yields 





:Tj|ir ^‘2- 1 ^ 




bii bii i^|ir 0 




(4i *** bi4 




^i! bis. (bi ir ) b>ai» 


ATp ^ 


ie 0 bj4 


] ~ 


b3t "biJe ba4 




^ _ 1 




C> Ui 0 _ t 



»• 

where 

A- 



A 



b., 




1 


0 


bi, 


b.2. 


bz3 


bL4 


bai 


bii 


0 


bj^ 


0 


b4^ 


b'43 


-I 



A 



Expanding and simplifying the notation results in: 

(I-iG) ~ ^ [jiz (c) ~-fi2 (B)J 



A' 



A' 



(X-17) 






^ \_hzi^) - kd^f\ 4 

A^ 



(A) — + bz 4 bij b +3 f bj2 — ^43 L34 

(B) — |_b'2.^ \mz -h b 3 i.-bi 4 .b 4-3 bi^.J 

( ^ ) ~ {^ ^ — t>i3 bu ~ bz4 kj 

(E ) - rbi.4 bsi b42. +-bii bsi — bit bai. - b4i bi4 bj. J 

1. B) =■ |_ ^31 biz, ~~ bii bai — b4^. bs4 bill 

fe) =[b,ib^» -b,, U^~b42, b^4 bu] 

Z\ — (■^)~ 4 i (B) A bj, (C )_J 



I 
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From ©(juations (1-16) and (1-17) it ia clear that the 
matrix elements are: 



/? = 

(1-18) 

(1-19) a,^ = 

(1-20) 5T,., = 



( 1 - 21 ) 




(/A) - -fzi ( 6) 





A' 






iii(c) - ■ 




(6) 




A' 








fe)- 


^21 


(F) 




A' 






4 


(&)- 


(Fj 



A' 



« 

I 

I 

J 
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APPENDIX II 



CALCULfiTION OF THE F-86 STABILITY DERIVITIVES AND 
DETERMINATION OF THE COEFFICIEOTS, AND f i j ,OF APPENDIX I 



A, Physical characteristics of the F-66, 

The following data has been obtained from Ref, <b . 
S = 268 ft? 

W = 15,500 lbs. (lightly loaded) 
b = 37.1 ft. 
c = 8.09 ft. 

I = 8625 slug-ft? 

I = 29400 slug-ft? 

I^x - 2700 slug-ft? 

1^2 - 35300 slug-ft? 

B. Assumed flight condition. 

Altitude = sea level 

Mach No, =; ,246 

density = ,002378 slug/ft^ 

Vp — 163 knots =.275 ft/sec 

h" ’*1= =o 

0(;l> — —3 = — 



The above data is for an F-86 proceeding down a three 
degree glide slope at 163 knots in a ’’clean** configuration. 

C. Lift and Drag Coefficients at Trim; Stability Dorivitives 



I. Cg 






577 
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2, From the drag polar using , 5?7 as the entering 

argument: Cn = . 04-4 V 

Therefore, Drag ~ 1152 lbs. 

ex', = 9.5° 

3, Stability Derivitives 

T^e following stability derivitives were obtained 
directly or indirectly from Eef. 6 . 



<k^ = .065 per degree = 3.72 per radian 

Ct), = 4^ = .561 per radian 

^ ~ < 7 d‘Ai — -.433 per radian 

P«r radian 



C 






.51 per radian 



=5.1 per radian 
C/^}^=-1.08 per radian 
=i «55 per radian 
= —.1775 per radian 
C^K = C/Vi^g .= -.945 per radian 






«.356 per radian 



^VC 

4. Coe fficients of longitudinal equations 

I htQ. CL-d cl I ~t~’ o nit- ^ Ct - ) 

Se, ^ z=. % = o 



7 



a 



• ^ . o/. ^ ^tr.M _ 

''-■j t ' "3 ^p. 



= — . OCO-5'4 



ACt- 

7f 



I 

I 



I 



^ ^ ■— xV " 



'5 





II 


— .(oO&>~I per 


■ dodiare. 




c. 


II 


± 








d. 




l^)cos>' 
-^J i*r/, • 


~ '3)^ 


<2. 


Ar'^v = 




II 


[~QVr.^] 




^2tJ< = 


)^p, 


■F,c,^J^0 




= 


— .7oQ -L 












/U>c. 








5* 




— ^ ^ / n &l0. 


~ 


'h. 00(pi3 — 












/K^C 


~k. 






— 


i€(- 


-^>j) = - 


* 

L. 






/^ 

( />i7p 


\ 

J =■ 


~ . oo6Z4 



J- = 



Cj, 

^ YTO ^(i 



(>Q. 2yr\ <PS<i 

* y'"”^ '' 

^ y4-^<lxXX-C^ / — O ^ 



o 



lb. i<^\/ 
1. Ka^ 



- O 






0/47 — , 

rad. 



.070 ~ 



yo. K<jo/ 



V. Kc 



n 



O. ^c 



r^ 



7 ""i 

ete! Cv,^ 



/yj 



. 3,09 -L 

-.‘oSl — 

■<^.74 
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D. Coefficients of Equations (1-9) through (1-12). 



11 



^12* 
^13 = 

bi4 = 
^11 = 
^21 = 
^22 “ 



= D/g - * 

= .6067 

^ 0 



.0311D -h .00654 



V = 

-K^v = .00085 
D- K2^= D +■ .709 
• _K^^ =-.00613 
‘’E4=-[i + Kzq]=-9S53 



'23 



^21 - 

^22 - 
bfj T ^ 



K ^ =+.00624 

zS-r 



'31 

b 7 o “ 



-V =0-00 



'32 

b^o = 



'33 
^34- 



- 

0 






.113D+ 3.09 



D - K = D 4- .567 

qq 



•31 



32 

^41 

^42 

^43 



44 



■41 



42 



= 0 

.Kq,^^.6.74 
= 0 

= D H d/at 

s d/at 

t= 0 
^ 0 
=i 0 



APPENDIX III 



ELEI.IENTS OF THE SEHVO AND TRANSDUCER MATRICES, T and S 
A. Servo Matrix Els ments, 

1. The relation between elevator deflection, X , and the 
signal to the elevator servo, ^ , is defined as 

A typical elevator servo system can be described as a 
second order system with natural frequency , UJn » and damping 
ratio, , The phase relation for such a system is given 




le ~ closed loop transfer function has a db 

versus log w characteristic as follows: 




It follows that if the elevator servo is designed such that 
the natural frequency, tOng , is larger than the expected 
range of system frequencies , Aui , or LA)«U)nethen equation 
(III-l) can be approximated by 

If the damping ratio is taken as .S^then the phase relation 
becomes : 

(III-3) (25 = 1 -4ti,n-' ~ -hxn~' n-'u) 



I 

I 



I 
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It may be shown that the phase relation of (III-3) is 
that phase relation for a simple first order system of the 
form : I _ / 

(III-4) 



In conclusion, if the damping ratio is chosen as ,5 , 

and the natural frequency of the servo is taken as 25 rad/sec^ 
-fer 

then'^the expected range of system frequencies the second 
order elevator servo system may be approximated by the first 
order system; — 



(III-5) o(p) 



£e 

4 



• 04 p 4- 1 



2, The relation between effective throttle perturbation or 
thrust, 4 » actual throttle perturbation, Sj- , is 

simply described by the following first order system: 



(III-6) 





where ^ = the time delay in obtaining the new 
engine EH.1 which is representative Of the 
thrust. 

The above relation is pointed out in Ref. 7 and is 
generally true for most turbojets. In addition, a typical 
time delay of two seconds , 'T'— 2 .0 , was obtained from this 
reference. Thus; 



(II 1-7) 




Zp-h 1 



B. Transducer Matrix Elements, . 

1. The transducer that relates actual indicated aircraft 
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velocity, /\J~p , to the measured aircraft indicated velocity 
is given by: 



Such a transfer function assumes no position errors in the 
airspeed indicator nor any dynamic errors in the instrument. 
The tirao*delay, 'T= .2 sec , is due to the pressure trans- 
mission delay in the pitot system. 

The basis of the foregoing statements is derived from 
a discussion of airspeed indicators and measuring devices 
contained in Ref, ^ . 

2. The transducer that relates actual glide path angle, 
Byj , to measured glide path angle, , is defined as tgg. 
It is assumed that the angle 0 and the angle of attack, ©<' 
are obtained without error from a gyro and an angle of 
attack indicator. These two quantities are added together 
to produce 0^, Therefore: 



(Ill-e) 




(III-9) 




± 






I 
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I 
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